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ABSTRACT 

A  surface  wind  (seabreeze),  thermally  generated  by  differential  sea  surface  temperature,  is  introduced 
to  Gill-Clarke’s  model  (1974)  through  wind  stress  for  investigating  the  effects  of  seabreeze  on  coastal  upwel- 
ling.  A  coupled  air-sea  system  is  treated  as  an  eigenvalue  problem.  The  solutions  show  that  the  thermally 
forced  local  winds  break  down  the  coastal  Kelvin  wave  into  three  parts:  small-scale  (L<100  km)  growing 
and  stationary  modes,  mesoscale  (100  km<£<200  km)  decaying  and  fast-moving  modes,  and  ‘large’-scalc 
(£>200  km)  coastal  Kelvin  modes.  The  consistency  of  the  length  scale  between  the  most  growing  mode 
predicted  by  this  model  and  the  observed  cold/warm  alternation  pattern  of  surface  water  near  the  Peruvian 
Coast  (around  15  °S)  implies  that  seabreeze  may  play  some  role  in  coastal  upwelling. 


I.  INTRODUCTION 

An  upwelling  study  JOINT  II  (Moody  et  al.,  1981 ;  Stuart,  1981)along  the  Peruvian  Coast 
(near  15’S)shows  that  sea  surface  temperature  (SST)  reveals  stationary  short  wave  pattern  in 
the  longshore  direction  (Fig.  1 ).  The  wave  length,  which  is  defined  as  the  distance  between  two 
adjacent  cold  (or  warm)  plume  centers,  is  around  55  km.  What  causes  the  formation  of  this 
small-scale  wave  pattern?  It  is  an  interesting  problem,  but  hasn’t  been  solved  yet. 

JOINT  II  data  also  show  that  the  area  and  temperature  of  the  plume  are  correlated  with 
the  wind  stress  (Stuart,  1981).  Up  to  now,  many  theories  treat  coastal  upwelling  as  a  wind- 
forced  problem,  i.e.,  the  surface  wind  is  taken  as  a  known  function.  In  fact,  surface  temperature 
gradient  generates  surface  winds,  which  in  turn  exert  an  extra  force  on  the  ocean  surface. 
Therefore,  coupling  the  air  with  ocean  is  important  in  investigating  coastal  upwelling.  In  this 
paper,  a  surface  wind  stress,  brought  on  by  thermally  forced  wind,  is  added  to  Gill-Clarke’s 
model  (1974).  The  model  starts  with  rest  ocean  and  atmosphere,  and  with  the  hori¬ 
zontal  homogeneity.  If  there  is  a  perturbation  initially  added  to  the  air-ocean  coupled  system, 
will  this  perturbation  be  growing  or  decaying?  be  oscillatory  or  nonoscillatory? 

The  subsequent  sections  are  tried  to  answer  those  questions,  to  prove  the  importance  of 
the  air-ocean  coupling,  and  to  predict  the  growth-rate  and  phase  speed  in  the  wavenumber 
spectrum.  The  length-scale  corresponding  to  the  maximum  growth-rate  and  minimum  phase 
speed  of  the  coupled  system  is  taken  as  the  general  appearing  length  scale  in  the  ocean  between 
the  two  adjacent  cold  (or  warm)  water  plumes. 
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Fig.  1.  Sea  surface  temperature  near  Peruvian  Coast  (after  Stuart,  1981). 


II.  AIR-OCEAN  COUPLED  MODEL 

In  this  study  the  air-sea  coupling  we  consider  is  shown  in  Fig.  2.  The  atmospheric  circula¬ 
tion,  thermally  generated  by  the  SST  gradient,  drives  the  water  current,  which  in  turn  changes 
the  SST  distribution.  The  salinity  effect  is  neglected  here.  Consider  the  ocean  of  uniform 
depth  H,  with  a  straight  coast  line.  The  right-handed  coordinate  system  is  chosen  that  the 
x-axis  is  in  the  alongshore  direction,  and  y-axis  is  in  the  cross-coastal  direction.  Furthermore, 
the  y-axis  is  nearly  coincident  with  1 5*S  latitude  (Fig.  3).  Therefore,  the  water  thermal  expan¬ 
sion  equation  is  written  by 

Pt->=rt->ci-«(r,-v7’v>):»  _  U) 

where  a  is  the  water  thermal  expansion  coefficient,  p\“\  TT0  are  the  basic  states 
of  water  density  and  water  temperature,  is  the  temperature  perturbation. 


A  state  of  rest  with  hydrostatic  equilibrium  and  horizon  tally-uni  form  stratification  is 
taken  as  a  basic  state,  i.e., 

dpvw)(z)!dz= (2) 
consider  the  theory  of  upwelling  for  the  case  where  only  small  perturbations  from  this  basic 
state  are  produced  (Gill  and  Clarke,  1974).  Let  the  perturbation  pressure  be  /»<“’,  the  den¬ 
sity  perturbation  be  a l“'),  and  the  velocity  components  be  («<“>,  a/1"’).  With 

the  Boussinesq  approximation,  these  quantities  satisfy  the  following  equations: 

( 3 ) 

fu"’'=-pr?'/ft'+Ymjifr'  a) 

0=  -  p\w>  -  golw)  (5) 

a)'1  +  o'/1  +  mV1’  =  0  (6) 

=  (7) 

where  the  subscripts  (/,  x,  y,  z)  mean  the  partial  derivatives  with  respect  to  those  independent 
variables,  and  /  is  time,  /  the  Coriolis  parameter  and  g  the  gravitational  acceleration. 
is  a  constant  which  is  a  representative  water  density.  The  vector  {X,  Y)  is  tangential  stress  act¬ 
ing  between  horizontal  planes.  According  to  Gill  and  Clarke  (1974)  we  have 


10,  z<~HmU, 

where  H„ix  is  the  depth  of  ocean  upper  mixed  layer.  Water  thermal  expansion 
equation  (1)  and  Boussinesq  approximation  give  a  relationship  between  water  density  per¬ 
turbation  a‘“)  and  water  temperature  perturbation 

Q,<”'  =  ~apk'")Tl’">  (9) 


X 


Fig.  2,  Air-sea  coupled  system.  Fig.  3.  The  coordinate  system. 


Author ’s  earlier  work  (1985,  1986ab,  1987ab)  shows  that  thermally  forced  surface  winds 
blow  across  isotherms  from  cold  to  warm  water  with  some  deflection  to  the  left  (right)  in 
the  Southern  (Northern)  Hemisphere  due  to  the  earth  rotation.  The  deflection  angle  <5  (Fig. 
4),  which  is  positive  (negative)  in  the  Southern  (Northern)  Hemisphere,  depends  on  the  atmos¬ 
pheric  stratification  (Chu,  1986b,  1987a).  |d|  is  generally  from  20°  to  30°.  The  more 

steady  the  atmosphere  is,  the  less  deflection  angle  S  should  be.  The  surface  wind  speed  is 
proportional  to  the  sea  surface  temperature  gradient.  Therefore,  the  thermally  forced  surface 
wind  is  computed  by 


y,»>=j*:roos<5  ~sin,5i  rar‘“7<3*-] 

Lsind  cosdJ  lSTm/dyX-*, 


(10) 


where K  is  a  proportionality  constant  which  represents  surface  wind  speed  thermally  driven  by 
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a  unity  surface  temperature  gradient.  According  to  Defant  (1950),  a  1°C/100  km  temperature 
gradient  can  generate  a  surface  wind  of  about  2  m/s.  Hence  we  take  (Chu,  1987b) 


K  =  ( 2  m/s)/(l‘C/lOO  km). 

(11) 

Using  following  bulk  formula  to  calculate  the  surface  wind  stress 

r<,‘>)=p.CcV'“\ 

(12) 

where  p„  is  air  density  at  surface,  and  CD  the  dimensional  (m/s)  air  drag  coefficient.  Sub- 

stituting  (10)  into  (12)  we  have 

rv\T'=pjCBK  (cosddT^/dx  -  sin 

(r*,0  ’  =  paCoK  (sin  8dTm/8x  +  cosddT^'/dy),^ . 

Gill  and  Clarke  (1974)  shows  that  there  exist  eigenfunctions  p„(z), 

(13) 

and  eigenvalues 

c„  such  that  u(u,>,  u'"1,  p"“\  and  Tm  can  be  expanded  as  follows: 

a*"1  =  2a*-"’1  (x,y,t)p„(z) 

(14a) 

(14b) 

p^/Pr  =  JlP^(=c,y,t)tn(z) 

(14c) 

(Hd) 

agTiw)/N !  (z)  =  'ZT^(x,  jr,f)  «M») » 

(I4e) 

where  JV  is  the  Brunt-Vaisala  frequency  defined  by 

N!(z)  =  -  ( g/pa)dp<,(z)/dz . 

(15) 

The  eigenfunctions  and  eigenvalue  for  barotropic  mode  are 

f  0  =  1 »  H)t  =  Z  +  H,  Cv=+/  gH  t 

(16) 

and  the  eigenfunctions  for  baroclinic  modes  are 
d'-ti>»/dz2+(N/c»)!ib!n=0, 

(17a) 

pn  —  dtiix/dz,  dfn/dz=  —  (N /ca)2ti)nr  n— 

and  satisfy  the  boundary  conditions  that  (Gill  and  Clarke,  1974) 

n  .(0)  —  0  J  «&.(--#)  =  0. 

(17b) 

Northern 
llem  i  spher  e 


Southern 
Hem i s  ph  e  r  e 

Fig.  4.  Thermally  forced  wind  and  surface  temperature  gradient. 


No.  3 


Instability  theory  of  air-sea  interaction 


«si 


=  E6.#.U>. 


The  forcing  terms  of  the  ocean,  Xt  and  Y2,  are  expanded  in  terms  of  the  eigenfunctions  ft„, 
e-g., 

r  x  j  i  r  t'/1  1 

[  (18) 
where  b„  («  =  0,  1,  2,...)  are  the  coefficients  of  the  expansion'. 

"«-{!:  -»»■«•  <•» 

Integrating  (19)  with  respect  to  r  from  -fitoflwe  find 

(20) 

From  the  temperature  expansion  (14e)  we  have 

Tlw>(x,y,0,t)=LN2(0)/ag^Ti,u>(x,y,t)^„  (0)  (21) 

and  from  the  top  boundary  conditions  for  eigenfunctions  **(’),  (16)  and  (17b),  we  obtain 

T<a‘>(x,y,0,t)=\:Nt(0)H/agin“\X,y,t).  (22) 

Substituting  water  density  perturbation  (9)  into  hydrostatic  equation  (5)  and  using  the  expan¬ 

sions  (14ce)  and  employing  the  relationship  between  eigenfunctions  p„  and  wn,  (17a),  we 
have 

T\r'=-p'?/cl.  (23) 

The  surface  water  temperature  perturbation  is  then  written  by 

rw(*,y,0,f)*-CAr,(0)/aflr,]pV',(*,y.O.  (24) 

Substituting  (24)  into  (13)  the  wind  stress  is 

<rlP= -ip,CnKN2(0)/agIl(cos}d  dp'pfdx-shnd  dpT/dy) 

-[paCnKN‘(0)/ag2l(smd  3p^/dx-rcos8  dp'^/dy)  .  (25) 

Therefore  the  basic  equations  for  the  n-th  mode  are 
du’-“)/8t  —  fv1,1  =  —  dp^jdx 

-y»(cos<5  dp'P'/dx -  s'ind  8p\f'jdy)  (26) 

—  dp^/dy  —  V»(sin5  dp\")/8x+  cosd  3p'“’fdy)  (27) 

dp^/dl  +  cl  (du^/dx  +  Sv\?'/dy)  =  0  (28) 

where  n=0,  1,  2,...,  and 

y»=V>  y^pV>C0KNHO)/(pim>Hmi2ag2) . 

Here  the  nondimensional  parameter  y  is  a  coupling  factor.  If  there  is  no  air-sea  coupling, 
y — 0,  the  system  is  reduced  to  the  ^Gill-Clarke  model. 


Ill,  INSTABILITY  CRITERION 

Eqs.  (26)-{28)  show  that  the  barotropic  mode  (n  -  0)  is  a  closed  system,  and  the  baroclinic 
modes  (n=  1,2,3,...)  are  forced  by  the  barotropic  mode.  If  under  some  circumstances  unstable 
barotropic  mode  is  generated,  then  all  the  baroclinic  modes  are  correspondingly  unstable. 
In  this  section  we  discuss  the  instability  criterion. 

Setting  n=0  in  (26)  (28),  substituting  the  differentiation  of  (25)  with  respect  to  y  from  the 
differentiation  of  (26)  with  respect  to  x,  and  using  the  continuity  equation  (28),  the  result  is 

dMC(i//)a4"’Va.v+/>r>/c3]  =  (y,sind//)v!f>r.  (29) 

Eliminating  from  (27)  and  (29)  then  gives 

dfdttl/P  8'-pl°’/8y2- pr/cl-rV'/P(sind  Vpf'/dxdy 
+  cos<5  a'j»r705>‘)]  =  -  (V„sin3//)V!JJr. 

Let  the  solution  have  the  following  form 


(30) 
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=$3<“',(:y)«x'p(0'f +  :£*) .  (31) 

Substituting  (31)  into  (30)  we  have  a  second  order  ordinaiy  differential  equation  for  ip'-w) 
d2q)('“)ldyl  +  LikV„  sin  6/(1  +  l\cos6)  jd<pt!U7c/y 

-  C/7  (1  -  Vooosd)  ]  (y.fe2  sin  8/ fa  +  1  /<*)*«  =  0 .  (32) 

Far  away  from  coast  perturbations  should  vanish,  i.e., 

as  y-*<x>.  (33) 

Right  at  the  coast  u‘“’  should  also  vanish,  i.  e., 

JvlT^du^/dt-dp^/dx  +  V'icosd  Sp^'/dx- sin6  8p^/3y)=  0, 

at  y  =  0.  (34) 

Eliminating  u'f  from  (27)  and  (34),  and  using  (31)  we  have  the  boundary  condition  for 

V<w) 

d<p(“‘ / d y  =  ik<pl“) (i  -*-Vi>cos6  —  opoSin  6/ f) /la If  J-70(sin6+cos6)3, 

at  ,y  =  0.  (3o) 

The  solution  of  (32)  satisfying  the  boundary  condition  (33)  is  written  by 

(pi'“)(y)  =  cexp(0y) ,  (36) 

where  0  is  a  root  with  negative  real  part  of  the  following  quadratic  equation 

0,  +  liky„sind/(l  +  V,<x)sd)10-fl(y,klamd/fa+l/ct)n\~V,ooBS)^.  (37) 

Substituting  (36)  into  the  boundary  condition  (35)  we  obtain 

/?  =  ife(l  +V(iCos6-  or„  sin  6//)  jla/f  +  y„  ( sin  6 -  cos 6)  ] .  (38) 

Eliminating  0  from  (37)  and  (38)  reaches  the  dispersion  relation 

IA\  -A,  +  pAJ  (a/}Y-lAt  +  2 A3  {A,  +  pAt)  ~r2A,A/]  (off)  '■ 

+  lAl+2A3At^Al(As  +  pA/)l(a/}')  +  AtA\=  0,  (39) 

where 


f.  =r/dk: 

=0.5y«sin6/(l -r '/doos6)  -y, 
J4,=l-!-y(,cos6J-0.5^sin6(sin6-cos6)/(l  -V„oos6) 
y4s=Vi,(sin6  +  cos6)  (40) 

A  4  ^  A  3  sin  6/  (l  +}'oCas6)  ~ 

A, s  -  0 . 2  5n  sin  ‘6/  ( 1  -  7tcos  6) 1 
^6  =  l/.(l+y<,cos6) . 

The  standard  values  of  parameters  representing  the  Peruvian  coastal  region  near  I5"S 
are  given  in  Table  1.  We  solve  the  characteristic?  equation  (39)  with  different  values  of 
k  and  6.  k  varies  from  10~3  m-’  to  KH  nr1,  and  6  from  20°  to  30°.  The  instability 
criterion  for  coastal  upwelling  in  the  air-sea  coupled  system  is  written  as 


!<0  decaying 

=  0  neutral  a~a„a3,a3,  (41) 

>0  growing 

where  a  is  the  root  of  cubic  equation  (39).  The  oscillation  criterion  for  coastal  upwelling 
is  given  by 


Im 


(a)  | 


=  0,  non-oscillatory 
oscillatory 


O'  —  O  1  f02>  O  3 


(42) 


The  cubic  equation  (39)  has  three  sets  of  roots:  first  one  contains  only  negative  real  roots 
for  all  k ,  which  represents  damping  non-oscillatory  modes.  The  other  two  sets  con¬ 
tain  complex  conjugates.  In  this  paper  we  are  interested  in  growing  modes,  so  that  we 
pay  less  attention  to  the  first  set  of  roots.  Furthermore,  without  the  air-sea  coupling 
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(v  —  0),  the  solutions  of  basic  system  (3)— (7)  should  be  coastal  trapped  Kelvin  wave  pro¬ 
pagating  in  the  alongshore  direction  with  the  coast  on  the  left  (right)  in  the  Southern 
(Northern)  Hemisphere  (Gill  and  Clarke,  1 974).  The  solutions  of  the  coupled  air-sea  system 
should  reduce  to  the  coastal  trapped  Kelvin  wave,  therefore,  we  should  choose  a  with 

Im  (<j)>0  (43) 

as  our  solutions. 

Table  1.  The  Standard  Values  of  Model  Parameters 

f  —  —0 .366  x  !0-4  5-1,  9  =  9.81  m  s~',  Cp  =  1.3  X  1(TS  ms”', 

tV2(0)  =  3,27x  10-1  s-2,  H  =  182.88  m  (10D  Fathoms) , 

p„  =  L.  29  kg  m-'’  9^.  =  103  kg  m-', 

0  =  2X10-*  K~',  11=2x10'  mV‘K-; 

Figs.  5  and  6  show  Re(cr)  and  Inri(o')  versus  the  negative  logarithmic  wave  number 
with  different  diffraction  angle  <3  (20°,  25°,  and  30°),  respectively.  They  indicate  fol¬ 
lowing  results: 

(a)  Thermally  forced  local  winds  produce  instability  in  the  coastal  upwelling  system. 
The  growth  rate  has  a  single  peak  in  the  wave-number  spectrum.  For  <3  =  20°  (very 
stable  atmospheric  stratification),  the  peak  is  located  nearly  at  k~  10=’-*,!i  m-1,  the 
corresponding  wave-length  is  L  =  2jr/fe  =  47  km.  For  <5  =  25°  (stable  atmospheric  strati¬ 
fication),  the  peak  is  located  at  fe=I0~‘-DTe  m_1,  the  corresponding  wavelength  is  L  — 
2jt/6  =  75  km.  For  <5=30°  (less  stable  atmospheric  stratification),  the  peak  is  located  at 
k  —  lO-1-156  m-',  the  corresponding  wave-length  is  L  —  2jr/fe=90  km.  It  shows  that 
the  wave-length  of  the  most  unstable  mode  enlarges  with  the  decrease  of  the  atmospheric 
stratification. 


Fig,  5.  Growth  rate  Re  (<r)  versus  wavenumber  K  for  three  different  deflection  angles : 
(a)  3=20°,  (b)  <5=25°,  and  (3)  3=30°, 
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Fig.  6.  Phase  speed  versus  wavenumber  K  for  three  different  deflection  angles: 
(a)  5=20°.  fl>)  5=25°,  (c)  5=30°. 


(b)  Fig.  5  indicates  that  according  to  the  growth  rate  distribution,  the  air-sea  coupled 
system  is  separated  into  three  parts:  small-scale  (£<100  km)  growing  modes,  meso- 
scale  (100  km<£<200  km)  decaying  modes,  and  large-scale  (£>200  km)  neutral  modes. 
However,  Fig.  6  shows  that  according  to  the  wave-speed  distribution,  the  coupled  system 
is  also  separated  into  three  parts:  small  scale  (£<(100  km)  stationary  modes,  mesoscale 
(100  km»  <200  km)  fast-moving  modes  (phase  velocity  is  greater  than  c0),  and  large- 
scale  (£(>200  km)  coastal  trapped  Kelvin  wave  (phase  velocity  is  c5).  Therefore,  the 
system  has  three  parts:  small-scale  stationary  and  growing  modes,  meso-scale  faster-moving 
and  decaying  modes,  and  large-scale  (£>200  km)  coastal  trapped  Kelvin  modes.  It  means 
that  the  coupling  effects  only  have  an  influence  on  the  small  and  meso-scale  modes,  not 
on  the  largescale  modes. 

(c)  Comparing  Figs.  5  and  6  with  Fig.  1,  which  shows  the  observed  stationary 
alternation  pattern  of  plumes  of  cold  and  warm  surface  Seawater  along  the  Peruvian  coast, 
we  find  that  length-scale  corresponding  to  the  maximum  growth-rate  (which  is  stationary) 
fits  the  length-scale  between  two  adjacent  cold  (or  warm)  plumes  of  the  surface  seawater. 
It  implies  that  thermally  forced  local  winds  may  play  some  role  in  coastal  upwelling. 

IV.  RESULTS 

This  coupled  air-sea  model  shows  that  the  thermally  forced  local  winds  break  down 
the  coastal  trapped  Kelvin  wave  into  three  parts:  small-scale  growing  and  stationary 
modes,  mesoscale  decaying  and  fast-moving  modes,  and  largescale  (£>200  km)  coastal 
trapped  Kelvin  modes.  The  alternation  pattern  of  cold  and  warm  plumes  of  surface 
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seawater  observed  by  Moody  et  al.  (1981)  and  Stuart  (1981)  can  be  predicted  by  this  model 
as  a  small-scale  growing  and  stationary  mode.  The  wave-length  corresponding  to  the 
maximum  growth-rate  also  fits  the  length-scale  of  the  two  adjacent  cold  (or  warm) 
plumes  very  well, 

In  this  model  we  neglect  the  nonlinearity  and  the  salinity  effect  on  the  water  density. 
We  should  include  these  effects  in  the  future  research. 

The  author  is  grateful  to  Professor  H.L.  Kuo  of  the  University  of  Chicago  for  invaluable  discussion  and 
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